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1. Introduction 

Numerical simulations in lattice field theory are only possible when both the lattice spacing 
and the volume are finite. If the thermodynamic limit is to be taken in theories with a mass gap, it 
is approached exponentially fast and usually the extrapolation causes no essential problems. Often, 
one is even interested in the universal continuum limit at a finite volume. It remains as the main 
problem to extrapolate the lattice spacing a to zero. Normally the use of some analytic formula 
for (at least) the asymptotic dependence of physical results on a — > is unavoidable 1 . Often this 
is based on the general framework given by Symanzik which for a universal dimensionless 
quantity P in many theories, including the a-model, amounts to setting 

a 2 

P{a)~P{0) + jjlco + c l \n{a/L)+c 2 \n 2 {a/L) + ..] +0(a 4 ), (1.1) 

where the coefficents c, have a finite continuum limit. In particular, in perturbation theory they con- 
tain only renormalized couplings, the bare one is eliminated before organizing the series in powers 
of In. For L we may take any finite renormalized length scale, an inverse mass for example. In finite 
volume calculations it is convenient to use the system extension. The derivation is based either on 
effective field theory ideas or on perturbation theory. There the power of the logarithm grows with 
the loop order. In applications in connection with Monte Carlo simulations one normally assumes 
that the function in the square bracket varies slowly over the range of accessible lattice spacings 
and it is often replaced by a constant. Then typically the form A + B{a/L) 2 is fitted to the data after 
discarding results for larger a/L until the fit is acceptable in terms of % 2 . Then the fit parameter A 
inherits an error 8 A from the data by error propagation, and in this way P(0) is estimated 2 . Clearly 
this procedure inflicts some systematic error, which may not always be negligible in high precision 
simulations as they are hopefully becoming more and more common in QCD in the future. We 
hence here take the occasion to investigate this issue in a simplified model where we have access to 
exact but still nontrivial nonperturbative results both in the continuum and at finite lattice spacing. 

In cluster simulations of the two-dimensional 0(3) nonlinear a-model [^, |l| over a rather 
large range of small a a behavior was empirically found for some quantities, which seemed hard to 
describe as a 2 but rather looked like a linear dependence. In [Q] we report additional simulations at 
N = 4 and N = 8 and a leading and subleading large N calculation. The main conclusion is that the 
often neglected In terms in ( |l.l| ) can and in fact do at large N mimic a behavior of the kind found 
in the simulations. We extend this investigation here by performing various extrapolations with the 
large N data and then assessing the systematic deviation from the known continuum limit. As a 
warning we shall demonstrate cases of significant discrepancy in combination with acceptable % 2 
and a plausible extrapolation formula. 

2. Step scaling function of the o -model at large N 

We start directly from the lattice formulation of the model where on a two-dimensional cubic 
1 See however ^ for an interesting attempt to avoid an explicit such step. 

2 Sometimes more conservative criteria are invoked at this point to exclude additional coarse data or add systematic 
errors. 
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lattice of size T x L the unit length A-component spin field s(x) with partition function 



Z 

is governed by the euclidean action 



fY\d N s8{s 2 -\)&- s ^ (2.1) 

J v 



^) = ^£(<V) 2 - (2-2) 

Here is the forward difference operator and y the rescaled coupling (held fixed in the large A" 
limit). 

As a physical observable that is finite in the asymptotically free (according to conventional 
wisdom) continuum limit we study the finite volume mass-gap [Q] M(L) that can be extracted on 
a torus with T — > °° at finite L from the decay of the fundamental spin correlation at zero spatial 
momentum. More precisely we study its scale dependence via the step scaling function (SSF) [§] 

L(u,N,a/L)=M(2L)2L for M{L)L = u (2.3) 

under factor-two rescalings of L with the continuum limit 

lim T,(u,N,a/L) = a(u,N). (2.4) 

a/L— >0 



The large A" evaluation of the path integral (2.1) is a standard saddle point calculation. More 
details on these steps can be found in [Hp. At leading order gaussian fluctuations of s(x) are con- 
trolled by the (euclidean) lattice Klein Gordon operator —A + mfc where the mass parameter mo is 
fixed dynamically in terms of y and a/L by the gap equation for T = °° 

H I IT! \ rrr^-i, ttj ' cb 2 = m 2 + p 2 , p = ^sm(a P /2). (2.5) 

y L p~^n 2( °(P)V 1 +a 2 (0 2 (p) /4 « 

The finite volume mass-gap is then closely related, sinh(a , M/2) = am /2 . 
If we expand the SSF as 

L(u,N,a/L) = L Q (u,a/L) + ^Li(u,a/L)+0(\/N 2 ), (2.6) 



the leading term Zo is obtained by twice using (2.5) - for some resolution L/a and the doubled 
value - and by combining the two to eliminate y. Evaluated numerically, this yields the exact 
SSF at finite cutoff. By an asymptotic expansion in powers of a using results of ^ we obtain 
the continuum value Ob( M ) an d the form of the leading artefacts. The limit is determined by the 
transcendental equation 

F(ob)=F(iO-ln2, F(z) = - + £{ * - - ) ■ (2.7) 



z „f i I \A 2 + (z/2n) 2 



n 



This continuum equation can be expanded in powers of u thus reproducing the large N limit of 
the known orders of perturbation theory of the j3 -function for this renormalized coupling [|5[ ](| 
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Figure 1: Correction at order l/N of the continuum step scaling function. 

|J. For large u the mass-gap becomes insensitive to the volume and Go(u) tends to the value 2u 
exponentially fast. The asymptotic lattice correction to On is rigorously known to be of the form 
( |L~I| ) with only cq,c\ non- vanishing. 

To obtain the 1 /N correction Ei we have computed the leading self-energy correction of the 
spin two-point function that shifts its pole by an order l/N term. The diagrams could only be 
evaluated [Q] numerically including a numerical integration. We managed however to compute 
them to many digits and could extract the continuum limit shown in Fig. [Ij The correction vanishes 
for both large and small u, since a is fixed uniformly in N in these limits. 

The investigation of the analytic behaviour of this correction is very hard. We know a repre- 
sentation of the continuum limit in terms of a numerical integral and a two-fold infinite sum which 
agrees with the lattice continuum extrapolation. However, it seems actually easier to perform the 
latter, as the sums also need an extrapolation. The small u expansion of the continuum form could 
be shown to agree with the l/N terms in the perturbative j3 -function up two three loop order. An 
extraction of the form of the asymptotic cutoff dependence has so far been beyond our capabilites. 
There is however strong numerical evidence [Q] that it is of the same type as in the leading term, 
i. e. just a 2 and a 2 lna. 

3. Sample continuum extrapolations of large N data 

We now use our 'exact' data for the nonperturbative values of Eo,£i to mimic their extrapo- 
lation based on lattices that could be simulated and with superimposed statistical noise. We test 
various more or less plausible fit functions. Our question is: when does this lead to significant 
biases beyond the error determined during the extrapolation? 

We consider two sequences of lattices for our virtual simulations 

if poor = {8,12,16,24}, Jgf rich = {10,12, 16,24,32,64,96, 128}. (3.1) 
For the lattices of either set we add artificial errors to £o 

Si = Io(m, Li/a)+wT]i, i = l,...,n, (3.2) 
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where {tj, } are gaussian random numbers of unit variance and w controls the size of the errors. We 
fit these numbers to a combination of n y independent fit-functions 

"/ 

Si « Fj = F(Li/a) = £ A a f a (Li/a), (3.3) 
a=l 

where we always take /i = 1 such that A i is the extrapolation. It is now a simple matter of linear 
algebra to show that on average (over the {fy}) we find 

= 1 + -TT- c££o(«, V«)^7 £ o(", V") ( 3 - 4 ) 

n — rif w L {n — nf)'rj 

with P the projector to the [n — n/) -dimensional space orthogonal to the one spanned by the 
f a {Li/a). We define that for a given set of functions a fit is accepted if Zmin/ ( n ~ n f) — 2 holds. 
Formula ( |3.4J ) allows now to determine a value wo such that this is true as long as w > wq. By 
similar algebra we find the mean propagated extrapolation error 



5A 1 =w^/(M- 1 ) 11 , M a p=Y,fa(L i /a)fp(L i /a) (3.5) 

i 

for the mean value of A\ obtained by solving the linear system 

YtMapAp =£/ (Va)Eo(«,V«)- (3-6) 

This in turn implies an error level, such that for w > w\ the systematic error is not significant 
in the sense \A\ — Gb(w) < SA\. In other words, for wo < w < w\ we typically find extrapolations 
that look acceptable but are intolerably biased. Of course, this is excluded, if wo > W\ holds. We 
now use our data for Eo an d Ei to determine the 'dangerous interval' [wo, w\] for a number of cases. 
We choose the value u = 2 where the cutoff effects, which are in general rather small at large N, 
are most noticeable for Eo- 0m - results are summarized in Fig. ||. The absolute errors w should be 
compared compared to Gb(2) = 3.2726 and Ci(2) = —0.0467. 



4. Conclusions 



We have sketched the computation of the step scaling function of the 0(/V) nonlinear a-model 
in the large ,/V limit including the subleading order. We focused on the approach to the continuum 
limit of these quantities which is reached with corrections of the form a 2 and a 2 In a. We used 
these data known beyond Monte Carlo precision to investigate various extrapolation methods, for 
instance in terms of pure powers. It turned out that for small enough errors there are cases possible 
with acceptable # 2 /d.g.f < 2 and an extrapolation bias beyond the statistical error of the extrapola- 
tion. As expected this becomes less likely if the assumed form is close to the true one and/or if the 
range of lattice spacings measured is large (except for the fit to a 15 ). As a very minimal precau- 
tion, one should try also forms with logarithms beside pure powers to probe systematic errors. The 
whole investigation may also be interpreted as a strong case for improvement programmes since 
flat data are easy to extrapolate. 
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Figure 2: Dangerous intervals for extrapolations for En (left) and Ei (right). The fit-functions employed are 



shown in the left of the figures, and solid and dashed lines refer to the sequences of lattices in (3. 1 ). 
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